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Abstract. 

We present the analysis of 200 high-resolution and 150 low-resolution spectra 
of the pulsating subdwarf close binary KPD 1930+2752 obtained with the Keck 
and the Calar Alto 2.2 m telescopes. Using metal-rich, line blanketed LTE 
model atmospheres the atmospheric parameters T c g = 35 500 ± 500 K, log g = 
5.56 ± 0.05 and log ^ = -1.48 ± 0.02 were derived. The radial velocity curve 
was measured and combined with all available published and unpublished radial 
velocity data in order to get results of high accuracy (K = 341 ± 1 kms -1 , 
P = 0.0950933 ± 0.0000015 d). For the first time we measured the projected 
rotational velocity. The preliminary result is v ro t smi = 97 ± 9 kms -1 (3c). 
Since the rotation of the sdB star is very likely tidally locked to the orbit, we 
can constrain the inclination of the system. Assuming the companion to be 
a white dwarf, the mass of the sdB is limited from 0.44 Mq to 0.63 Mq and 
the total mass range of the system is 1.3 Mq to 2.0 Mq. It is very likely that 
the total mass exceeds the Chandrasekhar limit. Hence KPD 1930+2752 is a 
candidate for a progenitor of a Type la supernova. According to the derived 
limits for the inclination angle, eclipses are likely to occur. 
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1. INTRODUCTION 

KPD 1930+2752 was identified as a subdwarf B star in the Kitt Peak-Downes 
survey (Downes 1986). The parameters which were derived from spectroscopy by 
model atmosphere fits (Saffcr and Liebert 1995) are consistent with the theoretical 
instability strip for pulsating sdB stars, which was predicted by Charpinet et al. 
(1996). After the first pulsating sdBs (EC 14026 stars) were discovered in 1997 
(Kilkenny et al. 1997; Koen et al. 1997; Stobie et al. 1997), Billeres et al. (2000) 
initiated a survey to search for these objects in the northern hemisphere. They 
selected KPD 1930+2752 from the list of Saffer and Liebert (1995) for their fast 
photometry program and detected multiperiodic variations with short periods and 
low amplitudes. In addition to 44 p-mode pulsations they found a strong variation 
at a much longer period of about 4100 s. This variation could be identified as 
an elipsoidal deformation of the sdB most likely caused by a massive companion. 
Billeres et al. predicted the period of the binary to be two times the period of the 
brightness variation (P = 8217.8 s = 0.095111 d). 

This was proven by Maxted et al. (2000), who measured a radial velocity curve 
of KPD 1930+2752 which could be fitted with the proper period. The radial ve- 
locity amplitude K = 349.3 ± 2.7 km s -1 combined with the assumption of the 
canonical mass for sdBs M s( jb = 0.5 M led to a lower limit for the mass of the 
system derived from the mass function. This lower limit M > 1.47M© exceeded 
the Chandrasekhar mass of 1.4 M Q . Because there was no sign of a companion 
in the spectra, it was concluded that the unseen object must be a white dwarf. 
Putting all this together Maxted et al. concluded that KPD 1930+2752 could be 
a good candidate for the progenitor of a Type la supernova. 

Type la supernovae (SNe la) are the most important standard candles for ex- 
tragalactic distance measurements and play an outstanding role in observational 
cosmology. The progenitors of SNe la as well as the dynamics of the explosions are 
still under debate. The thermonuclear explosion of a white dwarf turned out to be 
the most reasonable explanation. A white dwarf has to exceed the Chandrasekhar 
mass to explode in this way. In the single degenerate scenario the white dwarf ac- 
cretes matter from a nearby non-degenerate companion star. Possible progenitor 
candidates are supersoft X-ray sources. Double degenerate systems consist of two 
white dwarfs in a short period system, which merge within a Hubble time due to 
gravitational wave radiation (Livio 2000). KPD 1930+2752 has a sufficiently short 
period and may have a sufficiently high total mass. The sdB can be expected to 
evolve to a white dwarf before the merger. 

From the theoretical point of view Ergma et al. (2001) questioned the double de- 
generate scenario in the case of KPD 1930+2752. Their simulations based on the 
derived parameters of Maxted et al. suggested the formation of a single massive 
ONeMg white dwarf. 

The main drawback of all previous investigations is the lack of information on the 
inclination angle and the assumption made on the sdB mass. We drop the latter 
assumption and derive constraints on the inclination for the first time by means of 
an accurate measurement of the projected rotational velocity and surface gravity. 
Because the rotation of the sdB star is tidally locked to the orbit, we can derive 
sini as a function of the sdB mass from these two quantities. 
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2. OBSERVATIONS AND DATA REDUCTION 



With the 10 m Keck Telescope at the Mauna Kea Observatory two hundred 
high-resolution spectra were obtained by N. Przybilla in half a night in July 2004, 
using the High Resolution Echelle Spectrometer (HIRES, Vogt et al. 1994). The 
spectra covered a wavelength range of 4 200 A — 6 800 A with few small gaps at 
a resolution of 0.1 A and exposure times of 20 s each. The data were reduced 
using the ESO-MIDAS package. Bias and flatfield corrections were applied and 
a wavelength calibration was done. All spectra were corrected to the heliocentric 
frame of reference. 

Additional observations were obtained with the 2.2 m Telescope at the Calar Alto 
Observatory in July 2004. The Calar Alto Faint Object Spectrograph (CAFOS) 
was used to obtain 150 spectra covering a wavelength range of 3 600 A — 6 200 A 
with 5 A resolution and an exposure time of 180 s each. The data were reduced in 
analogy to the HIRES spectra. 



3. STELLAR PARAMETER DETERMINATION 



i 1 1 1 i 1 1 1 

He II 4686 

LTE log m = 



LTE log m = + 1 



The coaddcd CAFOS spec- 
trum had a very smooth contin- 
uum, which made it very suit- 
able for model atmosphere fits. 
Different LTE and NLTE model 
grids were used to fit 7 hydro- 
gen and 8 He I and Hell lines: 
A grid of metal-line blanketed 
LTE model atmospheres (Heber 
ct al. 2000) with solar and 
ten times solar mctalicity; a 
grid of partially line-blanketed 
NLTE model atmospheres with- 
out metal lines (Dreizler et al. 
1990) based on the ALI method 
(PR02); a grid calculated with a 
new version of PR02 (so called 
NGRT) that employs a tempera- 
ture correction scheme and uses 
more sophisticated model atoms; 
a small grid based on the LTE 
models with NLTE line forma- 
tion (Przybilla et al. 2006). The 
matching of the observed spectra was done with the FITPROF routine by the 
means of a x 2 fit (Napiwotzki 1999). The results of the fits with different grids 
were consistent (see Table 1). 
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Fig. 1. LTE models with logm = +1 solve 
the helium problem. 
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Table 1. 


Parameters of KPD 1930 


+2752 




model 


T c ff 


log g 


i n(Hc) 
l0 § „(H) 


LTE log to = 


35 549 K 


5.63 


-1.58 


LTE logm = +l 


35183 K 


5.56 


-1.48 


NLTE (PR02) 


35 619 K 


5.56 


-1.51 


NLTE (NGRT) 


35516 K 


5.53 


-1.45 


NLTE (If) 


35617 K 


5.64 


-1.53 


adopted 


35 500 ± 500 K 


5.56 ± 0.05 


-1.48 ± 0.02 



NLTE effects played a minor role. In the first instance it was not possible to fit 
the hydrogen and some helium lines (Hei 5876, Hen 5412 and in particular Hen 
4686) simultaneously. This so called helium problem already occurred during other 
analyses of pulsating subdwarfs (Heber et al. 2000; Edelmann 2003). The analysis 
of HST-UV spectra of three sdB stars with similar T eS as KPD 1930+2752 revealed 
supersolar abundances of the iron group elements (O'Toole & Heber 2005). Those 
stars also displayed the optical He ionisation problem. Using more appropriate 
metal-rich models (10 x solar metalicity) the problem could be remedied (see 
also Heber et al. 2006). The abundances of the iron group elements have not 
been measured for KPD 1930+2752. Because of the similarity of its atmospheric 
parameters to those of the stars studied by O'Toole and Heber (2005), we adopted 
high-metalicity models as well and, indeed, the fit improved (Fig. l).The derived 
temperature is about 2 000 K higher than in prior analyses (Billcrcs et al. 2000) 
irrespective of the choice of model atmosphere. With this KPD 1930+2752 is 
situated at the edge of the instability strip in the T e g — log g diagram (Charpinct 
et al. 1996). 

4. SPECTROSCOPIC ORBIT 

Due to the low signal to noise ratio of the HIRES spectra only the Ha and 
H/3 lines could be used for determining the radial velocity by x 2 cross correlation 
with a model spectrum at rest wavelength. To improve the accuracy, the resulting 
radial velocity curve was combined with all available radial velocity data of KPD 
1930+2752 (Maxted et al. 2000; Woolf et al. 2002; Orosz 2000 priv. comm.) 
covering a timespan of four years. A sine curve was fitted to these 2 900 data 
points using a x 2 minimizing method and the power spectrum was generated 
(FITRV and FITPOW routines by H. Drechscl). The sine curve fit is excellent and 
no period change could be detected over the whole timebase (Figs. 2). The orbital 
parameters were measured with unprecedented accuracy: 7(Ha) = 5 ± lkms -1 , 
K = 341 ± 1 kins" 1 , P = 0.0950933 ± 0.0000015 d. No radial velocity variations 
due to pulsations could be found. 

5. PROJECTED ROTATIONAL VELOCITY 

The main aim of the high-resolution time-series spectroscopy was to obtain a 
high-precision measurement of the projected rotational velocity. For this purpose 
the 200 HIRES spectra were shifted to rest wavelength and the median was cal- 
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Fig. 2. Radial velocity curve. The data points are from the new HIRES spectra, 
the solid line is from the fit of all data points over a timespan of four years. The 
residuals in the lower part show no signs of eccentricity. 
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Fig. 3. Reduced \ 2 as a function of v rot sini from fits to the HIRES and UVES 
data. The minimum of the two curves indicates the most probable value. The 
dashed horizontal lines mark the confidence levels of the errors. 
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culated in order to filter cosmics. The projected rotational velocity was measured 
using a model spectrum with the derived atmospheric parameters fixed and per- 
forming a x 2 fit- The lines Hell 4686, Hei 4922 and Hei 5016 were used for this 
measurement. Furthermore one single high resolution spectrum obtained with the 
ESO Very Large Telescope (VLT) at the Paranal Observatory and the Ultravio- 
let and Visual Spectrograph (UVES) was available. The rotational velocity was 
measured in the same way for the helium lines. The still preliminary results were 
consistent (see Fig. 3) and yield: v rot sini = 97 ± 9kms~ 1 (3a) 

6. MASS AND INCLINATION 

KPD 1930+2752 is obviously affected by the gravitional forces of the compan- 
ion, demonstrated by its elipsoidal deformation. Since the period of the photomet- 
ric variations are exactly half the period of the radial velocity curve the rotation 
of the sdB star is very likely tidally locked to the orbit. In addition to the mass 
function two more equations set additional constraints on the problem with only 
M s dB remaining as free parameter. 



With log g, P and ^ ro t sini measured and the fact that sini cannot exceed unity 
a lower limit for the mass of the sdB of 0.44 M Q is derived. As can be seen in Fig. 4 
the total mass of the system exceeds the Chandrasekhar limit for all reasonable 
assumptions of M s( jb ■ If the companion is a white dwarf, its mass has to be lower 
than the Chandrasekhar limit. This implies an upper limit M s( jb < 0.63 M Q and 
a possible total mass range of M s( jb+wd = 1-3 ± 0.1 — 2.0 ± 0.2 M Q (Fig. 5). 
Although a more massive companion cannot be ruled out completely, the object 
would probably be a bright X-ray source in this case. No such source is known at 
this coordinates. 

The inclination angle of the system (Fig. 4) in the lower mass range is very close 
to 90°. KPD 1930+2752 could be an eclipsing binary (cf. the very similar binary 
KPD 0422+5421, Orosz and Wade 1999), as some features in the light curve 
already indicated (Maxted et al. 2000). 

7. CONCLUSION 

The combination of an extensive set of 350 spectra and 2 700 radial velocity 
data points from prior observations made it possible to determine the parameters 
of KPD 1930+2752 with very high accuracy. The tidally locked rotation together 
with a white dwarf companion implies a mass range for the sdB without any 
evolutionary model considerations. KPD 1930+2752 is the first subdwarf, whose 
mass could be constrained in this way. The total mass and the merging time 
of the binary indicate that it is a very good candidate for a SN la progenitor. 
KPD 1930+2752 is one of only three known double degenerate systems, which 
fulfill these requirements (Napiwotzki et al. 2003; Karl 2004). While it fits into 
the DD scenario, KPD 1930+2752 might also evolve into a SN la via the SD 
scenario. Since time for merging due to gravitational wave radiation is of the 
same order as the EHB life time, Roche lobe overflow could occur well before the 




(1) 
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Fig. 5. Total mass of the binary as a function of the sdB mass. The dashed 
vertical line marks the point where the companion mass equals the Chandrasekhar 
mass. 



8 



S. Geier, U. Heber, N. Przybilla, R.-P. Kudritzki 



sdB becomes a white dwarf due to the shrinkage of the orbit. Recent calculations 
by Han and Podsiadlowski (priv. comm.) indicate that the mass transfer would 
be stable. Follow-up observations should be undertaken to measure an improved 
lightcurve and search for signs of an eclipse. If detected all system parameters will 
be determined. An alternative promising option is asteroseismology, which could 
tightly constrain the sdB mass. 

REFERENCES 

Billeres M., Fontaine G., Brassard P., Charpinet S., Liebert J., Saffer, R. A., 2000, 
ApJ, 530, 441 

Charpinet S., Fontaine G., Brassard P., Dorman B., 1996, ApJ, 471, L103 
Downes R. A., 1986, ApJS, 61, 569 

Dreizler S., Heber U., Werner K., Moehler S., dc Boer K. S., 1990, A%A, 235, 234 
Edclmann H., 2003, PhD-thesis, Friedrich Alexander Universitat Erlangcn-Niirnberg 
Eggleton P. P., 1982, ApJ, 268, 368 

Ergma E., Fedorova A. V., Yungelson L. R., 2001, A&A, 376, L9 
Han Z., Podsiadlowski P., Maxted P. F. L., Marsh T. R., Ivanova N., 2002, MN- 
RAS, 336, 449 

Han Z., Podsiadlowski P., Maxted P. F. L., Marsh T. R., Ivanova N., 2003, MN- 

RAS, 341, 669 
Heber U., 1986, A&A, 155, 33 

Heber U., Reid I. N., Werner K., 2000, A&A, 363, 198 

Heber U., et al., 2006, in Proc. of the Second Meeting on Hot Subdwarf Stars, 
Baltic Astronomy 

Karl C, 2004, PhD-thesis, Friedrich Alexander Universitat Erlangcn-Niirnberg 
Kilkenny D., Kocn C, O'Donoghue D., Stobie R. S., 1997, MNRAS, 285, 640 
Koen C, Kilkenny D., O'Donoghue D., Van Wyk F., Stobie R. S., 1997, MNRAS, 
285, 645 

Livio M., 2000, in Type la Supernovae: Theory and Cosmology, Cambridge Univ. 

Press, ed. Niemeyer J. C, Truran J. W., 33 
Maxted P. F. L., Marsh T. R., North R. C, 2000, MNRAS, 317, L41 
Napiwotzki R., 1999, A&A, 350, 101 

Napiwotzki R., Christlieb N., Drechscl H., et al., 2003, ESO Msngr, 112, 25 
Orosz J. A., Wade R. A., 1999, MNRAS, 310, 773 
O'Toolc S., Heber U., 2005, A&A, submitted 

Przybilla, N., Nieva M. F., Edelmann H., 2006, in Proc. of the Second Meeting 

on Hot Subdwarf Stars, Baltic Astronomy 
Saffer R. A., Liebert J., 1995, in Proc. 9th European Workshop on White Dwarfs, 

ed. Koester D., Werner K., Springer Vcrlag, 221 
Stobie R. S., Kawaler S. D., Kilkenny D., O'Donoghue D., Koen C, 1997, MNRAS, 

285, 651 

Vogt S., et al., 1994, Proc SPIE, 2198, 362 

Woolf V. M., Jeffery C. S., Pollacco D. L., 2002, MNRAS, 332, 34 



